arXiv:1505.07338vl [astro-ph.SR] 27 May 2015 


Binary Central Stars of Planetary Nebulae Discovered Through 
Photometric Variability III: The Central Star of Abell 65 

Todd C. Hillwig^’^i 

Todd.HillwigOvalpo.edu 
David J. Frew^’^ 

Melissa Louie^’® 

Orsola De Marco^ 

Howard E. Bond^’®’^^ 

David Jones^’® 

S. C. Schaubi’9’10 

ABSTRACT 

^Department of Physics and Astronomy, Valparaiso University, Valparaiso, IN 46383 USA and SARA 

^Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong 

^Department of Physics and Astronomy, Macquarie University, Sydney NSW 2109, Australia 

"^Department of Physics & Astronomy, Stony Brook University, Stony Brook, NY 11794 USA 

^Department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802 
USA 

®Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218 

^Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain 

®Departamento de Astrofisica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain 

^Southeastern Association for Research in Astronomy (SARA) NSF-REU Summer Intern 

^°Current address: Plasma Science & Fusion Center, Massachusetts Institute of Technology, Cambridge, 
MA 02139, USA 

^^Visiting astronomer, Cerro Tololo Inter-American Observatory, National Optical Astronomy Observa¬ 
tory, which are operated by the Association of Universities for Research in Astronomy, under contract with 
the National Science Foundation. 



2 


A growing number of close binary stars are being discovered among central 
stars of planetary nebulae. Recent and ongoing surveys are finding new systems 
and contributing to our knowledge of the evolution of close binary systems. The 
push to hnd more systems was largely based on early discoveries which suggested 
that 10 to 15% of all central stars are close binaries. One goal of this series 
of papers is conhrmation and classihcation of these systems as close binaries 
and determination of binary system parameters. Here we provide time-resolved 
multi-wavelength photometry of the central star of Abell 65 as well as further 
analysis of the nebula and discussion of possible binary-nebula connections. Our 
results for Abell 65 conhrm recent work showing that it has a close, cool binary 
companion, though several of our model parameters disagree with the recently 
published values. With our longer time baseline of photometric observations from 
1989-2009 we also provide a more precise orbital period of 1.0037577 days. 

Subject headings: binaries: close — planetary nebulae: individual (Abell 65) 


1. INTRODUCTION 


The complex shapes of planetary nebulae (PNe) provide a valuable and narrow evo¬ 
lutionary window through which we can study the life cycles of intermediate mass stars 
(typically 1-8 Mq). Binary interactions are known to play an important role but the extent 
to which such interaction s are needed to make a non-spherical PN is still unclear (for a 


review see 


De Marcoll2009[) . Studies have shown that approximately ten to twent y percent of 


all PNe harbor c/oge b inary stars with orbital periods typically less than 3 days flBondl 12000 


Miszalski et al.ll2009al). Such close binaries are the result of a common envelope (CE) inter¬ 
action (jPaczvnski 1976 '). The sample of known close binary central sta rs of planetary nebulae 
(CSPNe) more than doubled with the work of iMiszalski et al.l fl2009ar). However, stellar and 
syste r n parameters for t hese objects are very poorly known flDe Marco. Hillwig. fc Smith 
20081 : iHillwig et al.ll2010[ Pape rs I fc II). Unti l this si tuation is alleviated we cannot study 
this phenomenon statistically. iMiszalski et al.l fl2009bl) use this close binary sample to con¬ 
clude th at the nebular morph ologies of PNe with c lose binary central stars tend to be bipolar 
(see also Bond Sz Liviol 1990 1. Also, Hillwig (2010) provides a hrst look at the binary param¬ 
eters of the known systems and the relationship to similar systems which have no observed 
PN (either because they are more evolved, or because they did not form a visible PN at all). 
Additionally, there is increasing observational evidence for a class of intermediate-period 
CSPN binaries (with periods of r nonths to years). These binaries are discovered directly 
from radial velocity studies (e.g. IVan Winckel et al.ll2014j) . or otherwise inferred via evi- 
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dence of mass transfer from the PN’s progenitor sta r onto a companio n (e.g. iBond et al. 


2003; 

Miszalski et al. 

2012, 

2013b: 

Tvndall et al. 

2013) 

(2014 

)■ 


Planetary nebulae with post CE binaries are also extremely useful as probes of the CE 
interaction. Not only do they provide us with recent post CE binaries, where the period has 


not had time to evolve further, but they also show us the eject ed envelope itse 


they are extremely useful in understanding the CE efficiency (jPe Marco et ah 


also i n the timing of the various CE phases and the energetics of the ejecta flTocknell et ah 


As such. 


201 111, but 


2 OI 4 J) . For the latter purpose we need to know n ot only the details of the central binary, but 
also of the surrounding nebula (see for instance, I Jones et al.ll2014j) 


One of the goals of this series of papers is to further our understanding of the binary 
system parameters for known binary CSPNe. Connections between the binary parameters 
and the surrounding PN can then be explored, whether on a system-by-system basis, or 
through a statistical analysis of the larger sample. The former can be done as each model 
is produced. The latter depends on relatively few objects which currently have parameters 
derived from system modeling. We present here a study of the binary central star of Abell 
65, with a discussion of the surrounding PN and links between the two. 


2. BACKGROUND 


Abell 65 (PN G017 .3-21.9 l was first discovered bv ISharplessI fjl959l) and independently 
catalogued as a PN bv lAbelll (119661) . The 15th-mag central star was identified by Abell, 
who presented UBV photometry. Curiously, Abell did not note variability of the star, based 
on 4 separate observations, even though he did report variability for several other nuclei 
which was confirmed later. The variations of the Abell 65 nucleus were first discovered 
by H.E.B. in 1989. It soon became clear that the star was a close binary with an orbital 
period very close to 1 day, as was reported by Bond &: Liviol ( 1990 1. The ~1 day period 
probably explains why the large-amplitude variability was not detected by Abell: the orbital 
phase will be nearly the same at the same sidere al time every night during a short observing 
run! A lat er paper bv IWalsh fc WaltonI (1199611 erroneously classified the CS as an eclips¬ 
ing binary. IWalsh fc WaltonI (119961) also suggested that the CS was a cataclysmic variable 
(CV) s tar based on strong Ba l mer emission lines ob served in the spectrum. Improved light 
cur ves (IShimanskv et al.ll2009l: iLouie &: Hillwidl2010[l along with phase-resolved stellar spec¬ 
tra (jShimanskv et al.ll2009fl have since confirmed that the photometric variability is due to a 
strong irradiation effect, rather than eclipses, and that the irradiation effect is also the source 
of the strong emission lines observed in the optical spectrum. Those papers also confirmed 










































































4 


a period very close to one day. 


Physical values for the nebula such as distance , size, and age ar e given by a uum - 


ber of authors 


e.g. Walsh fc WaltonI Il996l: IPhillipsI l2005l: IPrewl I2OO8I: IPrew et ahl 120151) . 


Huckvale et al.l (120131) also provide a three-dimensional model of the nebula. They show 


a clearly 
tihed by 


jipolar nebula w ith both inner a nd outer structures ('the outer shell was hrst iden 


Huaetal 


ongmore 


19771). As not e d by iMiszalski et al.l (j2009bl based on an image from 


(jl9981) ). iHuckvale et al.l (120131) conhrm that while the bright nebular structure 


may resemble a toroidal waist, it is actually a set of bipolar lobes with fainter outer lobes 
(see their Figure 5). The inclination of the symmetry axis of these structures in their model 
are 68 ± 10° and 55 ± 10° for the inner and outer lobes, respectively. They also calculate 
kinematic ages for the inner and outer shells respectively, dependent on the distance to the 
nebula, of (8000 ± 3000) (d/l kpc) yr and 15, 000 ± 5000((i/l kpc) yr, whe re d is the dis- 
tance of the PN. D istances to the nebula range from d = 1.17 kpc from iFrewl (120081) to 
d = 1.67 k pc from (IStang he llini, Shaw, fc Villaverl 120081) . Hereafter we adopt the updated 
di stance ofiFrew et al.l ( 20151) . d = 1.4 kpc with a reported uncertainty of 20%. The values 
of iHuckvale et al.l (120131) and this distance give ages of the nebular shells of 11,200 yrs and 
21,000 yrs. 

To determine a redd ening value we star t ed with a we i ghted mean of the nebu l ar Balmer _ 

decrem ent e stimates fromlAcker et al.l (119921) , iKalerl (119831) , iPerinotto et al.l (119941) , IWalsh fc Walton 
( 1996 ), and Pollacco &: Bell (1997) to determine E{B — V) = 0.10. An independent value 
can be cal culated from a comp ariso n of the integrate d 1.4 GHz radio continuum and Ha 
fluxes from ICondon et al.l (119981) and iFrew et al.l (l2013l) . respectively; we do so to determine 
E{B — V) = 0.13 ±0.06. Also, the asymp t otic r eddening on this sight-line from the revised 
reddening maps of ISchlaflv fc Finkbeineii (j201ll) is E{B — V) = 0.09 ± 0.03. However, a 
strong IR signature in the PN from dust found in WISE and IRAS data suggests a poten¬ 
tially strong nebular contribution to reddening. This is conhrmed by lUE LWP and SWP 
spectra. The two spectra, taken at different times show clearly different fluxes, but neither 
match the slope of a h ot CS unless a r eddening va lue E{B — fo) > 0.18 is adopted. Using 
the reddening model of IValencic et ahl (12004120141) we hnd the best match for E{B — V) = 

0.20 and adopt this as the reddening value for our modeling. 

In Figure m we show the reddening corrected lUE spect ra, U, R, and I fluxes fr om this 
work, and XMAS'S* ( Skrutskie et ah 2006 ). WISE, and IRAS ( Taiitsu fc Tamura 1998 ) fluxes. 

Taiitsu fc Tamura ( 19981 ) use corrected IRAS fluxes to determine a dust temperature of 65 
K. Their correction assumes extended dust emission over the entire optical nebula, so they 
scale the emission up from the size of the aperture to the full size of the nebula. Here we 
simply use the reported values, two of which are upper limits, as shown in the hgure. 
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Using the ephemeris in §4.1 we find that the SWP spectrum was obtained at phase 
0.90, or close to minimum light (which occurs at orbital phase 0 = 0). The LWP spectrum, 
obtained at 0 = 0.78, has greater fiux than the SWP spectrum, suggesting that part of the 
heated hemisphere of the companion is hot enough to affect the near-UV spectrum. We find 
that this requires a temperature of at least 20,000 K on some portion of the companion. 
Our final reddening value quoted above was found using the SWP spectrum along with our 
minimum light V and I values. The remaining data from 2MASS, WISE, and IRAS were not 
obtained at minimum light so we do not expect them to match the included model, but they 
are shown to demonstrate the IR-excesses due to the companion and dust in the nebula. In 
particular, the WISE and IRAS bands use large apertures and so include both thermal dust 
emission and fine-structure line emission, as well as light from the central binary. Indeed, in 
the WISES and WISE4 band images, the nebula overwhelms the star in brightness. In the 
three IRAS bands, the magnitude difference between the observations and our stellar/dust 
model are 0.82, 0.99, and 1.09 magnitudes in the J, H, and K bands respectively. The IRAS 
data were obtained at orbital phase 0.77 based on our ephemeris given below. For that 
orbital phase, these values are slightly higher than we would expect based on the irradiation 
effect alone. Likewise, the WISEl and WISE2 band values are 1.37 and 1.55 magnitudes 
above our minimum-light model, higher than expected from irradiation. However, these 
values are within a reasonable range when nebular contamination is also considered. We 
label these data points in Figured] as upper limits to remind the reader that due to nebular 
contamination and not being obtained at minimum light, they are upper limits for the model 
shown in the plot. The WISES and WISE4 bands are used, along with the IRAS fluxes, to 
estimate the dust component of the system. 


Without higher resolution in the infrared it is difficult to determine the source of the 
dust emission. However, given that our best fitting reddening value is significantly higher 
than values from a number of other methods as described above, it seems likely that there 
may be an additional compact dust source around the CS that contributes to the reddening, 
such as the dust disks reported for a number of PNe flBilikova et al.ll2012l. e.g.). In Figure 


[T]we show a two cool dust components: a 115 K blackbody to fit the WISES and WISE4 
bands, and 40 K blackbo dy fit to the IRAS 6 0 /rm a nd 100 pm bands. The cooler component 
roughly matches with the iTaiitsu &: Tamural (119981) dust emission from the nebula, while the 
warmer component is similar to dust disks seen in these systems, which average around 150 
K. In both cases, it would appear that there is an additional reddening contribution by dust 
in the nebula. 


In §4 we further discuss several of the literature values for the Abell 65 nebula in the 
context of our binary star modeling. 
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Wavelength (/xm) 


Fig. 1.— The spectral energy distribution plot for the CS of Abell 65. The lUE spectra, 
V, R, and I magnitudes [solid circles) at minimum light, 2MASS J, H, and K magnitudes 
[open squares), WISE fluxes [open triangles), and IRAS fluxes [open circles). All fluxes have 
been de-redenned assuming E[B — V) = 0.20. Also shown are model blackbody curves for 
our resulting CS, companion, and dust component, as well as the combined curve. These are 
shown without any irradiation contribution, thus should approximately represent the system 
at minimum light. 
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OBSERVATIONS AND REDUCTIONS 


The photometric data consist of images obtained over twenty years at fonr different 
telescopes in the V, R, and I hlters. The Cerro Tololo Interamerican Observatory (CTIO) 
0.9 meter and 1.3 meter telescopes, the SARA 0.9 meter telescope at Kitt Peak National 
Observatory, and the Perth Observatory 24-inch telescope were all nsed to collec t data. The 
data inclnde the original discovery observations described in iBond fc Livid (1199011 and follow- 
np observations by H. E. Bond on the CTIO 0.9 meter telescope. These data were collected in 
1989 September, 1990 Jnne, 1991 September, and 1992 May (U, R, and I hlters). The SARA 
0.9 meter telescop^ was used to make observations from the northern hemisphere at multiple 
epochs in 2007 and 2008 {R hlter). The CTIO 0.9 meter and Perth 24-inch telescopes were 
used to observe the CS of Abell 65 on three consecutive nights in 2008 August, providing 
coverage of almost 70% of the orbit each day (U, R, and I hlters). The CTIO 1.3 meter 
provided observations in queue observing mode from August to October 2008 and June and 
July of 2009 {R hlter). The Perth 24-inch telescope also provided observations in queue 
observing mode from 2008 June-August (/ hlter). 

All images were bias subtracted and hat helded (and dark subtracted, when necessary). 
The IRAF/DAOPHOT package was used to perform aperture photometry on the images. 
The 2007-2008 obse rvations were the n jointly analyzed using incomplete inhomogeneous en¬ 
semble photometry (iHonevcuttl Il992h . The ensemble photometry used all of the stars in a 
template image to perform diherential photometry. The large number of stars of various col¬ 
ors produces an ehective color correction for slight diherences in system throughput between 
the four diherent telescope, hlter, and camera combinations. The success of this method 
was conhrmed in the hnal light curve by the seamless overlap of the 2007-2008 data from 
the four telescopes in the phase-folded light curves (Figure [2]). The 1989-1992 data were 
analyzed by single star diherential photometry with two check stars. The earlier data were 
then shifted to the same magnitude zero-point as the more recent data (see §4.1 below). 
Some of these earlier data do not align well with the 2007-2008 data, especially in the V 
hlter. These diherences could be a result of diherent CCD and hlter combinations, or could 
be intrinsic to the system (see §4.2 for more discussion). The two data sets were combined 
to provide a more precise orbital ephemeris. 

The resulting light curves are in an instrumental magnitude system. However, the 
photometric data were then shifted to apparen t magnitudes us i ng th e V = 15.80 and I = 
15.39 magnitudes from the HST observations of ICiardullo et al.l (119991) . The hnal ephemeris 


^Based on observations obtained with the SARA Observatory 0.9 m telescope at Kitt Peak, which is 
owned and operated by the Southeastern Association for Research in Astronomy (saraobservatory.org). 












for the binary system was used to determine the phase of the HST images (0.81 and 0.79 
for V and I, respectively) and our photometric data were shifted to those values at the 
determined phases based on our hnal binary model light curves. Given the high precision of 
our ephemeris, we adopt a magnitude zero-point uncertainty of 0.06 mag based on the small 
uncertainty in orbital phase and the slope of the light curve at that phase. To determine the 
shift necessary for our R hlter data, we use a log £t between the V and / de-reddened fluxes 
to extrapolate an R magnitude value for Abell 65. The resulting zero-point is thus slightly 
less well constrained, but should still be reliable to within 0.1 magnitudes. The result gives 
apparent magnitudes at minimum light oi V = 16.04, R = 15.80, and I = 15.66. These 
minimum brightness values provide the best observable measure of the combined brightness 
of the two stars in the system with a minimum irradiation contribution. They are used with 
the lUE SWP spectrum to arrive at a value for interstellar reddening, as shown in Figured! 


The spectrum of the CS of Abel l 65 was obtained o n 2008 July 11 with the UVES high 
resolution spectrograph on the VLT flDekker et al.ll2000l) . The exposure time was 1200 s with 
a spectral resolution of R ~60,000 and covering the range of 4800-6800 A. The spectrum 
was obtained at an airmass of 2.285 with the slit at the parallactic angle of 107°. The seeing 
was 1.38 arc seconds. The spectrum was reduced with the standard ESO UVES pipeline. 
Details of the spectrum are discussed in §4.3. 


4. ANALYSIS 


4.1. Orbital Ephemeris 


Because we were unable to analyze the early epoch data using the same method as the 
more recent data (see §3 above) the two data sets had magnitude scales with different zero 
points. To correct this a preliminary model of t he binary sys t em was produced using o nly 
the recent data with the Wilson-Devinney code fjWilsonlll990l : I Wilson fc Devinnevlll97ll ) in 
the V, R, and I hlters. The resulting model was then £t to the early epoch data, allowing 
only the magnitude zero-point to change. The derived zero-point offset was applied to the 
early epoch data to shift them to the same zero point as the more recent data. Having all 
of the data on the same zero point scale allowed us to calculate the orbital ephemeris to our 
highest possible precision. To do so we used the Wilson-Devinney code to £t the photometric 
data with only the time at phase 0 = 0, Tq, and the orbital period. Port, as variables. The 
resulting binary system ephemeris is then 


T = 2454301.3244(38) + 1.0037577(10) x 0 
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where 0 = 0 is taken to be at inferior conjunction of the cool secondary star (what would 
amount to primary eclipse if this were an eclipsing system). This phase also corresponds 
to minimum light. The values in parentheses reflect the uncertainties in the last two digits 
of each number. The complete phase-folded V-, R-, and /-band light curves are shown in 
Figure [21 where each has also been shifted to the apparent magnitude scale as described 
above. 


4.2. Binary System Parameters 


Shimanskv et ah! (12009!) use their spectra and P-band light curve to calculate a model 
for the binary system confirming that it is indeed an irradiated binary with a hot central 
star and cool companion. However, we believe that with our light curves in the V, R, and I 
bands along with some additional information described below, we are able to produce a more 
accurate model of the system. Because of the large number of parameters, binary system 
modeling often requires that a number of those parameters be estimated independently and 
used as fixed va l ues in the modeling in order to reduce the degeneracy of the solution. 
Shimanskv et ah! fl2009l) use the size an d expans i on rat e of the PN to arrive at a nebular age 
and then use the evolutionary tracks of iBlockerl (11995!) to derive a radius of the central star. 
The primary difference in our approach to modeling the system occurs at this point. With 
our determined apparent magnitudes at minimum light, the published distances to Abell 65, 
and interstellar absorption values we determine the absolute magnitudes for the combined 
light from the binary CS of Abell 65. We then constrain our models to match these values. 

Using minimum light in the V band {V = 16.04), the smallest distance value in the 
literature (d = 1.17 kpc, iFrew! 12008!) and our adopted reddening (giving Ay = 0.62), we 
hud = 5.09 for the CS in Abell 65. Using the largest value of distance, d = 1.67 kpc 
(IStanghellini. Shaw, fc Villaver! [2008I) we hnd an alternative r ninimum ab s olute magnitude 
for the CS of My = 4.31. And with our adopted distance from iFrew et ah! (!2015[) of d = 1.4 
kpc. My = 4.69. 


The binary parameters of [Shimanskv et ah! (1200911 result in an absolute magnitude. 
My = 3.2, well outside our range of 4.55 < My < 5.33. It appears that the major cause of 
the discrepancy comes from their calculation of the age of the nebul a (and thus the radius o f 
the CS). They quo te a radius for the exten ded halo of Abell 65 from !Walsh fc Walton! (1199611 
of 2 pc. However, !Walsh fc Walton! (jl996!) use a distance of 1.5 kpc to a rrive at a diameter 
for the nebula of 2 pc rather than a radius. The [Shimanskv et ah! (12009!) resulting age then 
is a factor of two too large. Using an as s umed expansion velocity and a larger distance than 
our adopted value of 1.4 kpc from [Frew! (12008!) leads to additional discrepancies in age. The 
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result is that their v alue for the age o f the nebula of 130,000 years is considerably different 
than the result from Huckvale et ah fj2013 ) of 21,000 years and 11,200 years for the outer 
and inner nebular shells, respectively (for qur ad opted distance). Even with the difference 
in age it is not clear how IShimanskv et ahl (120091) arrive at their primary radiu s of 0.12 Rcd. 
Extrap olating between r nasses of 0.546 Mq and 0.565 Mq from the models of ISchonberner 
(1l983[) used by iBlockerl (119951) one could find a value of the mass for which the CS radius 
would be approximately 0.12 Rq, but the basis for their claim that it is the most likely radius 
is not made clear. However, their final model mass is 0.56 Mq for which an age of 130,000 
years places the CS much closer to the WD cooling track with a radius approximately a 
factor of hve smaller. It is largely due to the sensitive relationship between CS mass, radius, 
and age that we have added the absolute magnitude as a constraint on our models. However, 
as discussed below we do use theoretical evolutionary tracks to help guide our modeling. 


Using the Wils o n-Dev inney code with our photometric data and the radial velocities 
of IShimanskv et ahl (120091 ) we find a best-fit solution for the system parameters. In order 
to limit the possible parameter space we use the absolute magnitude co nstraint described 


above along with the temperature of the cool companion star found from IShimanskv et ah 


(120091) from their spectra, T 2 = 4950 K. The best-fit values from our modeling are given 
in Table [Hand over-plotted with the data in Figure [2l Because the uncertainties in htting 
the data in binaries such as this are dominated by systematic effects, the uncertainties given 
for each parameter are not statistical errors for that individual parameter, but represent the 
range over which each parameter can vary in correlation with the other parameters. The 
resulting uncertainties are larger than those found by holding all other parameters constant 
while varying only one, but provide a more physically significant range for each value. 


In order to properly fit the light curves we departed from using, for the secondary star, 
the internally calculated limb-darkening values in the Wilson-Devinney code. We found that 
lower values were necessary to reproduce the light curve shapes and amplitudes. The model 
shown in Figure [2] uses the internally calculated limb-darkening values for both the R and 
I bands, but the V band value has been reduced slightly. The change in limb darkening 
required to fit the observed light curves for any combination of parameter values in the 
ranges given in Table [T] is within 10% for all three bands. 

We note several discrepancies in the photometry visible in Figure [2l Both the V and 
R data on the increasing brightness side fall slightly above the model, though the effect is 
most evident in V. The asymmetry is real and while most of the data in that portion of the 
light curve are from a single night, there are a few data points from three different nights, 
suggesting that the discrepancy is not instrumental or due to atmospheric conditions. We hnd 
that we cannot match the discrepancy with our modeling. The source of the asymmetry could 
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potentially be due to non-synchronous rotation of the secondary and some of the irradiated 
hemisphere rotating around toward the back hemisphere. The greater effect with decreasing 
wavelength, the effect is not noticeable in the I data, would suggest an increase due to higher 
temperatures, as would be expected from asynchronous rotation. The orientation is also in 
the correct sense for prograde rotation of the secondary at a rate greater than the orbital 
period. There is no corresponding dip in the light curve below our model on the decreasing 
brightness side of the light curve, as might be expected if synchronous rotation is moving the 
irradiated portion of the atmosphere. However, an asymmetry in that portion of the curve 
would only occur if the irradiation time scale is on the order of the rotation time scale. If 
the irradiation timescale is much shorter, then no such asymmetry would be expected. We 
cannot at this point determine whether non-synchronous rotation is the culprit, but it would 
seem to produce the observed features. Another possible explanation would be a bright spot 
on the cool companion, though a consistent bright spot at the same longitude on the star 
would seem unlikely. 

We also find slight discrepancies in the peak of each light curve. In R and I the data were 
taken on different nights with the Perth telescope and the small brightness change appears 
to be inherent to the system. Finally, the V data on the decreasing brightness branch from 
the early CTIO data (obtained in the 1990’s) lie clearly below the model while the early R 
data fall above the model. The I data from the 1990’s line up with the more recent data. 
It is possible that a difference in filter response is responsible. Or there may be real changes 
in the amplitude over time. 


The light curve amplitudes from our model are AV = 0.882 mag, AR = 0.912 mag, and 
AI = 0.914 ni ag. The relative amplitude difference follows the relationship we presented in 
Paper I ( 2008 ). 


The resulting absolute magnitude of the system from our model values is My = 4.70 


which agrees well with t he value found above (M y = 4.69) using the distance of iFrew et ah 


(120151) . As discovered bv IShimanskv et al.l (120091) . a hot and bright CS is required to £t the 
data, and we find that we are able to fit the data for relatively large ranges of temperature 
and radius of both the CS and companion. However, the degeneracy in the models is reduced 
dramatically by using an assumed system brightness. 


In additi on to our modele d abso lute magnitude, several other system parameters differ 
from those of IShimanskv et al.l (120091) . We find a companion with a lower mass, M 2 = 0.22± 
0.04Mq and radius, R 2 = O.41±O.O5i?0 compared to their values of M 2 = 0.33±0.06Mq and 
R 2 = 0.59 ± 0.08i?Q, respectively. As with many other cool secondary companions in close 
binary CSPNe, we hnd that the secondary in Abell 65 has a higher temperature and radius 
than expected for a main sequence star of similar mass. Thus the secondary is over-luminous 
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Fig. 2.— The V, R, and I phase-folded light curves of Abell 65 for the period given in the 
ephemeris (Equation 14. ip . 
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by about a factor of 20-30. Our M 2 value suggests a companion spectral type M4V while the 
temperature would be equivalent to a K2V star and the radius is appropriate for an M2-3V 
spectral type. So both the temperature and radius are large compared to those expected for 
a main sequence star based on our modeled mass. 

The temperature of the CS of 110,000 K is larger than, but consistent with th e Zanstra 
temperatures we calculate using the integrated Ha flux from iFrew et al.l fj2013|) and the 
dereddened Hell/Ha ratio averaged from the references given in Section 2. The results are 
Tz H = 49 ± 7 kK and He = 93 ± 12 kK (the large discrepancy between the hydrogen and 
ionized helium Zans tra temperatures i ndicat es the nebula is optically thin). It is also slightly 
higher than that of IShimanskv et al.l (120091) . but the two agree to within the uncertainties. 
However, we hnd a CS radius more than a factor of 2 smaller than theirs. As a consistency 
check we use our resu l ting M cs, logfTd,^) , and log(Lc 5 '/L 0 ) with the evolutionary 

tracks of ISchonberneii (1198311 and iBlockerl (1l995l) . The results are highly consistent with 
those tracks and lead to a CS age of about 22,000 years which, within the uncertainties in 
m ass, radius, temperatur e, and distance, is in good agreement with the nebular age estimated 
bv iHuckvale et al.l (1201311 of 21, 000 ± 7, 000 years (for the outer nebular shell, scaled to our 
distance). The value of 11, 200 ± 4, 200 years from iHuckvale et al.l (120131) for the inner shell 
is younger than the age we hnd from the evolutionary tracks. This may be signihcant in that 
we would expect a CS tha t underwent CE ev olution to evolve much faster than a non-CE star 
(cf. Ricker fc Taam 2012 : Passv et ah 2012). So if the CS in Abell 65 began its post-AGB 
evolution at a time consistent with ejection of the inner nebular shell, then we see a much 
shorter ac tual age than the pre dicted evolutionary age assuming a single star. Because the 
models of iHuckvale et al.l (l2013l) are based on spatio-kinematic model we consider their ages 
to be the most reliable. Even for the largest published distance for Abell 65 the age of the 
inner shell is still younger than the evolutionary age of the CS. We do note however, that 
small deviations in the CS mass will produce a large change in the evolutionary age, making 
this comparison of limited quantitative use given our uncertainty in Mcs- In addition, the 
starting points for the evolutionary tracks are somewhat arbitrary and the spatio-kinematic 
models assume no acceleration or deceleration of the expanding nebula. So we again state 
that using this age comparison only serves as a consistency check for our results, rather than 
an absolute comparison. 


Our inclination, i = 61 ± 5°, is lower than that of IShimanskv et al.l (120091) [i = 
68.0 ± 2.0°), as a result of the different stellar parameters. Our value though is i n bet - 
ter agreement with the nebular inclination for the inner shell from iHuckvale et al.l (120131 ). 
i = 55±10°, and al s o agre es with their value for the inclination of the outer shell, i = 68 ±10°. 
As IHuckvale et al.l (120131) state, the connection between close binary CSPNe and their host 
nebula is becoming clearer. For the handful of systems that have had both nebular and 
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binary in clinations determi ned, the tw o values are within the uncertain ties in every case: 


Abell 41 (Jones et ah 

2010). Abell 63 (Mitchell et ah 2007). NGG 6337 (Garcfa-Diaz et ah 

2 OO 9 I: Hillwie et al.lbolO). HaTr 4 jTvndall et ah 20121. also Hillwig 2015a in prep). NGG 6778 

(Miszalski et ah 

2011; 

Guerrero & Mirandal 20121). Henize 2-428 (ISantander-Garcfa et ah 

2 OI 5 I). and Sp 1 

(Jones et ah 2 OI 2 I also Hillwig 2015b in prep). 


4.3. The Spectrum of the Central Star of Abell 65 


In Figure [3] we show a high resolution nebular subtracted spectrum of the CS of Abell 65 
(see §3 for details) in the regions around H/5 and Ha. The spectrum was obtained at orbital 
phase 0 = 0.255 based on our ephemeris. What immediately stands out are the very broad 
H I lines with centr al absorption. These lines are reminiscent of those seen in the CS of HFG 1 
(lExter et al.ll2005l ) with very broad double peaked H I lines. The line shape is reminiscent 
of the double-peaked H I profiles of accretion disks in cataclysmic variables (CVs), thus 
the occasional classification of HFG 1 and Abell 65 as GVs. Ho wever, this b road double- 


peaked line prohle is well described for irradiated atmospheres by iBarman et al.l (120041) with 
the central absorption due to non-LTE effects in the upper atmosphere of the irradiated 
companion star (see their Figure 7) . The consistent change in equivalent width of these lines 
through the orbit (see Figure 4 in IShimanskv et al.l 120091) further conhrms that the origin 
is not an accretion disk, for which the line strength would only ch ange appreciably in the 
event of an eclipse. Furthermore, the minimum light spectrum from IShimanskv et al.l (120091) 
shows the hot star to be H-normal, giving an inferred spectral type of 0(H). The strong 
irradiation effect is also evidence of the non-GV nature of Abell 65 since for a GV the system 
brightness would be dominated by the disk, thus swamping the irradiation effect and greatly 
reducing the observed variability amplitude. 

We measured the radial velocity of the central absorption of the H I lines in the VLT 
spectrum and the narrow G II emission lines at 6578 and 6583 A, both of which originate from 
the secondary star. The central absorption in H I could be effected by incomplete or over 
subtraction of the nebulosity, but we find radial velocity values consistent with the narrow 
G II lines, suggesting that any nebular sub traction effect is s mall. The resulting velocity 
point is plotted in Figure 0] with the data of IShimanskv et al.l (120091) and our best-fit model 
parameters. The VLT radial velocity point agrees well with both the previous data and our 
model. 
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Orbital Phase (2455301.3244 + 1.003757?) 


Fig. 4.— The phase-folded ra dial velocity plot for the CS of Abell 65. Shown are the data 
from IShimanskv et all (120091 ) {solid circles), our single data point {open circle) from the VLT 
spectrum at phase 0.255, and our model £t {solid lines) which includes our predicted curve 
for the hot CS. 
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CONCLUSION 


Shimanskv et al.l (120091) have shown conclusively that the CS of Abell 65 is an irra¬ 
diated close binary CSPN. Due to the large parameter space available for binary system 
modeling, determining a hnal parameter set typically requires more information than radial 
velocity and light curve s. With our three-color photometry we hnd that the parameters of 
Shimanskv et ahl (120091 1 do fit well. However, we believe that our use of luminosity as an 
initial fitting requirement results in a more accurate and more consistent parameter set than 
their initial requirement of CS radius based on nebular age and CS mass. Our results do 
not change the overall identification of the system as a hot central star and cool, irradiated 
companion. However, we believe that our final set of physical parameters more accurately 
reflects the make-up of the binary system. And we believe that what may seem relatively 
small differences in some of these parameters are important for future interpretations of the 
increasing body of close binary CSPNe. 

Obtaining precise physical parameters for these systems can allow us to explore several 
interesting areas. First, we find that in cases where both the nebular inclination and binary 
inclination are known, the values show a very strong correlation. As the number of systems 
with known inclinations increases we will be able to say more about the likelihood that close 
binary interactions dominate the shaping mechanism in the resulting PN. Second, as we 
increase the statistical significance of the close binary CSPN sample we can begin to explore 
the birthrates of cataclysmic variables (CVs) and type la supernova progenitors of various 
types. Based on our results, the binary CS of Abell 65 will reach the semi-detached stage at 
a period of approximately 0.29 days and beco me a CV in ap proximately 9 x 10® years due 
to magnetic braking, based on calculations in IWarnerl (119951) . The time to a possible single 
or double degenerate Type la SN is longer than a Hubble time, as is expected for systems 
with periods of more than about half a day. 

As is the case for the large majority of post-CE CSs, this binary CS in Abell 65 has a 


period smaller th an three days. Common envelope simulations (e.g., iRicker fc TaamI 12012 


Passv et al.1120121) predict a broader range of post-CE binary separations than is observed. 
It is this type of comparison that has alerted us to the fact that our model of the CE phase 
is far from complete. The irradiation properties of the companions in post-CE binaries are 
also going to serve as constraints of the CE interaction. Here too, CSs of PNe are unique, 
because the bright and hot primaries result in the brightest irradiation effects. One parameter 
that may be constrained is the amount of accreted angular momentum on the companion. 
Irradiation theory should be able to tell us whether the companion is synchronized or not, 
and hence give us information on its rotation properties. These properties connect to both 
its magnetic field, which is currently held responsible for bright and hard X-ray emission 
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(e.g., 


Montez et al.l I 2 OIOI: iKastner et al.l 1201211 and the amount of accret ion that may take 


place before or during the CE (e.g., iMiszalski. Boffin, fc Corradill2013al) . Finally, as with 


most other secondary stars in close binary CSPNe, the secondary star here is both heated 
and expanded, thus over-luminous, compared to a main sequence star of similar mass. This 
characteristic of close binary CSPNe is likely a result of the CE phase, though it is not clear 
why the stars have not returned to thermal equilibrium. Increasing the statistical sample 
of over-luminous secondaries will hopefully lead to a better understanding of the CE and 
post-CE phase effects on the companion star. 


This material is based upon work supported by the National Science Foundation under 
Grant No. AST-1109683. Any opinions, findings, and conclusions or recommendations 
expressed in this material are those of the author (s) and do not necessarily reflect the views 
of the National Science Foundation. This publication makes use of data products from the 
Wide-held Infrared Survey Explorer, which is a joint project of the University of California, 
Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded 
by the National Aeronautics and Space Administration. This publication makes use of data 
products from the Two Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Center/California Institute of 
Technology, funded by the National Aeronautics and Space Administration and the National 
Science Foundation. 

Based on observations made with ESO Telescopes at the La Silla Paranal Observatory 
under programme ID 081.D-0857(A). 
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Table 1. Abell 65 Binary System Best-fit Model 


Parameter 

Value 

q 

0.40 ±0.10 

Tcs (X10^ K) 

110 ±10 

T2 (x103 K) 

4.95 ±1.0 

Mcs (Mq) 

0.56 ±0.04 

M2 (Mo) 

0.22 ±0.04 

Res (Re) 

0.056 ±0.008 

R 2 {Rq) 

0.41 ±0.05 

Igfi-CS (cgs) 

6.7 ±0.2 

lgfi'2 (cgs) 

4.6 ±0.2 

^(°) 

61 ±5 

a (Rq) 

3.90 ±0.10 

log(-^cs/ Lq) 

2.61 

log(T2/Lo) 

-1.05 






